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Theoretical equations are developed for typical decompositions of polymers including those in 
which the volatilization does not follow a simple "reaction order" and those made up of a composite 
of several reactions of differing energies of activation. The effects of order, activation energy, heating 
rate and temperature dependence upon the calculated thermograms is illustrated. The literature on 
thermogravimetric kinetics is critically reviewed and coalesced into a logical and coherent development 
stressing the interrelation of methods and employing a consistent system of notation. As a result, 
a number of improved methods and new methods for the analysis of kinetic data applicable to the 
complex systems mentioned above are developed. It is concluded that methods involving a variable 
rate of heating or involving several thermogravimetric traces at different rates of heating are capable 
of establishing the uniqueness of kinetic parameters. A new method of determining initial parameters 
from rate-conversion data is developed. A novel concept is employed of programming reaction variables 
(in this case, the heating rate) in a manner which greatly simplifies the mathematics of the kinetic 
system and which shows promise of a wide range of applicability in the area of rate processes. 
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1 . Introduction 

Thermogravimetric analysis is used widely as a method to investigate the thermal decomposi- 
tion of polymers and to assess their relative thermal stabilities [1, 2,3,4]. l Also, considerable 
attention has been directed recently toward the exploitation of thermogravimetric data for the 
determination of kinetic parameters. A number of these methods will be discussed later in this 
paper. 

Many of the methods of kinetic analysis which have been proposed are based on the hypothesis 
that, from a single thermogravimetric trace, meaningful values may be obtained for parameters 
such as activation energy, preexponential factor and reaction order. 

Thus, many of these methods make two assumptions, viz, these parameters are useful in 
characterizing a particular polymer degradation, and that the thermogram for each particular set 
of these parameters is unique. 

Therefore, in this paper, we will test the validity of these assumptions by setting up several 
idealized typical cases of polymer degradation kinetics, observing how the structures of the cal- 
culated thermograms are affected by changes in order, activation energy, heating rate and tempera- 
ture dependence, and determining by means of a critical evaluation of both existing and new 
methods if there are any general treatments of these data that will permit the determination of 
parameters that may be useful in the interpretation of degradation mechanisms. 

2. Theory 

We shall assume for the present that the isothermal rate of conversion, dC/dt, is a linear 
function of a single temperature-dependent rate constant, k, and some temperature-independent 
function of the conversion, C, i.e., 

f=*/(Q. (1) 



1 Figures in brackets indicate the literature references at the end of this paper. 
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C, the conversion (degree of completion or advancement, extent of reaction), is defined here as 
the conversion with respect to initial material. Thus C=l — (WIWo), where Wo is the initial 
weight and W is the weight at any time. Therefore, the residual fraction (1 — Q = W/Wo and the 
rate of conversion, dC/dt = — l/Wo (dW/dt). 

Some investigators prefer to define an instantaneous rate of conversion such that dC /dt 
= \\W {dWjdt) and thence C is proportional to the logarithm of the residual fraction. However, 
unless perhaps dealing specifically with first order reaction kinetics, there seems to be little 
special appeal for this latter definition of conversion. 

When the polymer does not completely volatilize at T or £— > oo 9 or if the volatilization takes 
place in steps, the conversion may for convenience, or depending upon one's insight into mech- 
anism, be defined differently based on the total weight loss between two successive horizontal 
portions of an integral thermogram. Equation (1) excludes composite cases where simultaneous 
or successive reactions involve several temperature dependent rate constants. Several of these 
cases will be considered later. 

At a constant rate of heating, jQ = dT/dt, and assuming k independent of C and/(C) independent 
of T, the variables in eq (1) may be separated and integrated to obtain: 

C c dC 1 C T 

where F(C) is the integrated function of conversion and 4> represents the temperature integral. 
In analogy to simple cases in homogeneous reaction kinetics, we express the conversion func- 
tion by 

/(0=(n+l)(l-0» (3) 

where n is defined as the order or reaction, n + 1 is the normalization factor for the isothermal 

cases such that (n+ 1) I (dC/kdt) dC=l. n is equal in value to n and included only in the cal- 

Jo 

culated curves so that they will exhibit maximum thermogravimetric rates at approximately the 
same temperature. 

Substituting eq (3) into eq (1) and integrating, one obtains 

-=-fa or-$ (n *1) 



(l-n)(n + l) 

1 



(4) 
In (1-0 = - At or -O (n=l) 



(n+1) 

(-fa, if T= const.; -$, if 7V const.). 

The normalized isothermal curves from eq (4) for rate versus conversion are plotted in figure 
1. Some of these curves indeed do approximate ones obtained in polymer degradation studies. 
For example, if the depolymerization is initiated at chain ends and the zip length of the depropaga- 
tion reaction is much shorter than the polymer chain length, as is the case in high molecular weight 
polymethyl methacrylate, then a large portion of the reaction follows zero order kinetics [5]. If 
the zip length is much larger than the polymer chain length, first order kinetics results [5]. Other 
cases of degradation kinetics, such as is found for branched polyolefins, may be fitted over a range of 
conversion by higher order curves [5]. 

However, in a large number of polymer decompositions, the isothermal rate of conversion goes 
through a maximum value. Such behavior describes a wide spectrum of polymers [5]. 
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FIGURE 1. Effect of order on normalized isothermal rate 
versus degree of conversion. 
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Random depolymerization, where most of the early scissions result in fragments too large to 
distill off, produces an isothermal rate curve with a maximum at 26 percent conversion. 
For random depolymerization [6], 



l-C=(l-a) 



1 + a 



(N-L)(L-l) 



N 



(5) 



where IV — 1 is the initial number of carbon atoms in the chain skeleton, L is the number of carbon 
atoms in the smallest chain that does not evaporate, and a, the fraction of bonds broken. 
Differentiating eq (5) with respect to kt, 



dC =(L-w-c)-(i- a y- {N - L * L - l) 



(6) 



dkt 



N 



and eliminating a from eqs (5) and (6) for the simplest case, N> L, L = 2, result in 

f(C) = 6(C^-C) 

for the normalized isothermal equation. In formulating the above approximation (Case A), the 
term for the initial isothermal rate, [L(L+ 1)]//V has been assumed equal to zero since in practice 
N > L. Hence the expression goes to zero at C = Q and is not applicable to initial situations. 
The solution of eq (5) for a involves an Lth order polynomial so that it is not practicable to obtain 
f{C) in analytic form for large L. Also, maximum rates are never obtained at greater than 26 
percent conversion from random depolymerization. Therefore, an empirical cubic equation was 
devised to fit several cases found in isothermal degradation kinetics. 



For the equation 



dC 
dkt 



= N(C 3 + a 1 C 2 +a 2 C + a: i ) 



(7) 
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the following boundary conditions: 




\dkt) c =o> \dkt) max ' 


(dC\ \dkt)c=o> 
\dkt)c-i ' /dC\ 7 
\dkt)c=o 


are satisfied when 


' 


12[(y-l)(l-2co)-^] 


_2a, 3 -(y-l)(l-3o> 2 ) 

' fl| /. . 1 Wl OX 9 '» 



w^txt)^ — loo) — yj — [y — i) low — oo) — ij [y—i){i — zo)) — co' 

= qj[(y-l)(2-3aj)-aj 3 ] co 2 (co-l) 2 

tt2 (y-l)(l-2aj)-a> 2 ' anda3 (y -i)(i_2 w )-a) 2 ' 

For Case B, we assume the maximum rate at 50 percent conversion and a maximum rate three 
times the initial rate, (a; = 0.5, y = 3) and obtain, 



*>—!!(■ 



+c 2 -|c-^ 



For Case C, (co = 0.25, y = 3), eq (7) has a minimum as well as a maximum in the range C = 
to 1. However, Case C was fitted satisfactorily by reduction of the order of eq (7) which gave after 
normalization, 



w=-fi{ cm+c -l cl,2 -\)- 



The normalized isothermal rates for Cases A, B, and C are plotted against conversion in 
figure 2. 

Curves A and C closely resemble theoretical curves for pure random depolymerization or 
for depolymerizations of moderate zip length involving considerable molecular transfer. Case B 
is representative of depolymerizations in which the initiation occurs at chain ends, the zip length 
is moderate and either a slight amount of transfer or random initiation take place. Such curves 
are characterized by a maximum in the 50 percent conversion region [7]. 

The equation for Case A may be substituted into eq (1) and integrated and if the solutions of 
the cubic in the equations for Cases B and C are expressed as partial fractions, substituted into 




FIGURE 2. Normalized isothermal rate versus degree of 
conversion, Cases A , B, and C. 
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eq (1) and integrated, one obtains: 

for Case A, 

3 \n(l-C l ' 2 ) = -kt or-<S>, 

for Case B, 

| {ln(l-C) +0.6547 In (C + 1.8660) -1.6547 In (C+0.1340)- 1.6219} =-A* or-O, 

and for Case C, 

y|r {In (1-C 1 / 2 )- 1.2217 In (C^+l. 8660) -0.2217 In (C 1 ' 2 + 0.1340) + 0.5451} =- to or -4> 

(— kt, if T= constant; — 4>, if T t^ constant). 
The Arrhenius equation, 

i=^-« (8) 

is almost universally assumed for the temperature dependence of A:, where A, the "preexponential 
factor," is usually assumed to be independent of temperature, E is the energy of activation, and 
R, the gas constant = 1.987 cal mob 1 fK)" 1 (1 cal = 4.1840 J). 
Substituting eq (8) into eq (1), one obtains 

§- fi §-A+*"m (9) 

which upon integration becomes 

where x = { — E/RT). It is assumed that T is low enough for the lower limit to be negligible. 
If /(C) is given by eq (3) then eq (9) becomes, 

^T = P ^=Ae^ E l RT \l-C) n (11) 

at dl 



which, upon integration becomes, 



(l-CV^-l AE , , 



l-n $R 



(n^l) (12) 

AF 

\n(l-C)=^-p(x) 

(» = 1). 

p{x), which includes the exponential integral, has been tabulated for limited ranges [8-11]. 
There are several series expansions and a semiempirical approximation for p(x). These are given 
below as they are utilized in many of the kinetic methods which will be discussed later. 
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2.1. Asymptotic Expansion [12] 
px / 91 V 4 ! 



x ^ — 10; number of terms ^\x\ . (13a) 



2.2. Expansion in Series of Bernoulli Numbers [12] 

/% e * ( aaaa^aoc , 0.998710 , 1.98487646 , 4.9482092 , 11.7850792 
pyx) = — — . 0000035 + — 1 = 1 h - 



/y \ y /y^ yO /V"^ 

-A^ \ ,A/ A/ A/ -A/ 



+ 2 M ^ + 2 M 4|l^ + 9 : 52^411 ±0 _ 35xl0 _ 



2.3. Schlomilch Expansion [13] 

e -y / 1 2 4 14 

P(y) = Mr+l)V~rT2 + (y + 2)(y+3)~(y+2)---(y+4) + ( r +2)---(y+5)~' 

(y=-*^15). (13c) 

2.4. Doyle's Approximation [14] 
log p(*) = - 2.315 + 0.457* 

(-20^*^-60). (13d) 

Doyle [15] has found that two and three term approximations for eq (13c) are more accurate 
than the respective approximation for eq (13a). Equation (13b) is almost equivalent to eq (13a) 
for the first several terms. 

At the limits, — 20 ^ x ^= — 60, eq (13d) is accurate within ±3 percent [11]. 

An extensive table of log p(x) for various x's was calculated from eq (13b) as this equation 
gives quite accurate values even for small x. This table was used for the calculation of the ther- 
mograms presented subsequently. 

If the preexponential factor, A, is linearly dependent on temperature, i.e., A = A t T 

* = AiW-*/« (14) 



then, 



yhere, 



•--^g+Hl? *}--*£** <*> 



^ e * ( nnnnMK 00010290 , 1.9487646 , 
pi(x) = — 1—0.0000035 1 - 1- . . . ). (16) 



A table of log pi(x) values was composed to investigate the effect of temperature dependence. 
Vallet [9] has derived expressions for A = AzP 12 , k = A 2 7 n ' 2 e (-£! ""' ) and recursion formulas 
for calculating the above three cases of temperature dependence. 
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3. Calculated Thermogravimetric Curves 

Integral thermograms were obtained by calculating <t> and T at various x and substituting 
4> into eq (17). 



(l-0 = {(n + l)(^-l)0-l-l} T=H (n¥> 1) 
(1-C) = e ( - 2(I,) (n = l). 



(17) 



Figure 3 contains curves for the residual fraction, 1 — C, versus temperature for various orders 
from to 4 for the case A 1/3 = 10 16 /°K, E = 60,000 cal/mol. 

The differential thermograms corresponding to the integral curves in figure 3 were obtained 
by substitution of values of <$> at various T into 



dC = (n + l)A { (n -n )(yi - !)$+!} n-i e (-MV 
dl B 



(18) 



These thermogravimetric rate versus temperature curves are contained in figure 4. 

Zero and negative orders exhibit ever accelerating rates. Curves for orders between zero 
and unity go through a maximum rate and, in the idealized case, reach completion at a finite 
temperature. For curves with n ^ 1, increasing order results in a more gradual slope and a more 
gradual asymptotic approach to the abscissa. 

Integral thermograms for the "maximum" Cases A, B, and C were obtained for the same 
values of A//3 and E by calculating 4> at various values of C from the integrated equations for the 
three cases. The corresponding x and T values were obtained from <J>. Residual fraction versus 
temperature curves for these cases are contained in figure 5. The early negative-order character 
of these curves manifests itself by a more precipitous slope than for the zero order curve at the 
same activation energy. 




Figure 3. Effect of order on residual fraction versus 

temperature. 

n = 0, 0.5, 1.0, 2.0, 3.0, and 4.0; 

Alfi" 10 ,6 /°K; £ = 60,000 cal/mol. 
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FIGURE 4. Effect of order on thermogravimetric rate 
versus temperature. 

« = 0, 0.5, 1.0.2.0, 3.0, and 4.0; 
/<//3=10 1 «/°K; £=60,000 cal/mol. 
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Figure 5. Residual fraction versus temperature for Cases 

A, B, and C. 

/*/£=10 16 /°K; £=60,000 cal/mol. 



FIGURE 6. Thermo gravimetric rate versus temperature 
for Cases A, B, and C. 

Alf3= 10 16 /°K; £ = 60,000 cal/mol. 



The corresponding differential thermograms for Cases A, B, and C were obtained by direct 
substitution of C and the corresponding T values from the integral calculations into eq (9). 

The thermogravimetric rate is plotted against temperature for these cases and the same values 
of A//3 and E in figure 6. The first order curve is included for comparison. 

The minor detail of the differences among Cases A, B, and C is not significant as these curves 
are based on asymptotic or semiempirical equations. The slope is greatest in Case B where the 
maximum isothermal rate is at 50 percent conversion. As y (maximum rate/initial rate) approaches 
unity, the curves would more nearly resemble zero order with a tail. 



3.1. Thermogravimetric Rate Versus Conversion Curves 

The previous figures (figs. 3-6) of residual fraction or thermogravimetric rate versus tempera- 
ture have been the traditional methods of representing thermogravimetric data. It is surprising 
that plots of the rate as a function of degree of conversion (dC/dT versus C) have not been utilized 
as similar presentation has been found to be quite informative in isothermal studies, e.g., figure 1. 

Figure 7 exhibits the thermogravimetric rate as a function of conversion for zero through fourth 
order and Cases A, B, and C for A/ (3 = 10 16 /°K and E = 60,000 cal/mol. 

The zero order curve is almost linear while higher order curves go through a maximum rate 
as C— >1. They approach the zero order curve as C^0. Curves B and C also approach the 
zero order curve asymptotically at low conversion but deviate in the direction of negative order. 
Curve A is anomalous at C^0 as was mentioned previously. 

A better insight into the behavior of the curves in figure 7 may be obtained by differentiating 
eq (9) with respect to conversion to obtain 



d 2 C E f'(C)dC 
dCdT Rr f(Q dT 



(19) 



where /'(C) is the derivative of /(C) with respect to C. If /(C) is given by eq (3), then similar dif- 
ferentiation of eq (11) gives 



dC 



(PC E 
dCdT RT 2 1-Cdf 
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(20) 



FIGURE 7. Thermogravimetric rate versus degree of 



dC 
dT 



conversion. 

n = 0, 0.5, 1.0,2.0.3.0, and 4.0 

Cases A, B, and C 

A/p= lO'VK; E =60,000 cal/mol. 




Thus the slope of the zero order curve will be given by E/RT 2 . At ordinary activation energies, 
the absolute temperature is a slowly changing variable over the reaction temperature range so the 
slope will be essentially a constant. The thermogravimetric rate is quite small at low conversion 
so all other cases approach the zero order case as C— >0, subject to the condition, f(C = 0) / 0. 

Thus, one may obtain a qualitative understanding of the early kinetics from the limiting 
characteristics of dCjdt versus C. If the curve is initially concave, a negative order character, 
i.e., an initially increasing isothermal rate as in Cases B and C is indicated. Initial convexity 
suggests positive order behavior, while approximate linearity connotes kinetics not far removed 
from zero order. Indeed such a plot may be of great practicable applicability, e.g., to test whether 
an apparent maximum rate in an isothermal curve is due to an initial temperature lag or an actual 
initial increase in rate. 

3.2. Effects of Activation Energy 

The effect of change in energy of activation upon the character and shape of the thermograms 
was investigated by considering three cases of first order kinetics— I, £ = 60, 000 cal/mol, II, 
£" = 40,000 cal/mol, and III, £ = 80,000 cal/mol. The parameters A//3 were adjusted so that the 
maximum slopes occur at the same temperature in each case. The integral curves are plotted in 
figure 8 and the corresponding differential curves in figure 9. 

The maximum slope increases linearly with increasing energy of activation since, from eq 
(20), at the maximum 



dC\ 
dT) n 



£(1-C) n 



(21) 



Other relationships at the maximum will be discussed in the section on differential methods. 
3.3. Effect of a Temperature Dependence of the Preexponential Factor 

A slight temperature dependence of the preexponential factor is sometimes noted in isothermal 
studies. Indeed, according to collision theory, A is proportional to T 1/2 for a bimolecular gas 
phase reaction and, according to transition state theory, the preexponential factor contains tem- 
perature to the first power. 



495 





FIGURE 8. Effects of activation energy and preexponen- 
tial temperature dependence on residual fraction versus 
temperature. 
n=l 

I £=60,000 cal/mol;/<//3 = 10 i rK 

II £=40,000 cal/mol; /4//3=9.675x 10 9 /°K 

III £=80,000 cal/mol; ^//8=9.318x 10»/°K 

IV £=60,000 cal/mol; A/0= 1.374 X 10 13 /°K; A{T\ 



Figure 9. Effects of activation energy and preexponential 
temperature dependence on thermogravimetric rate 
versus temperature. 
n=l 

I £=60,000 cal/mol; Alfi= 10 18 /°K 

II £=40,000 cal/mol; y<//3=9.675x 10»/°K 

III £=80,000 cal/mol; A/ft =9.318 x 10"f K 

IV £=60,000 cal/mol; /<//3=1.374x 10 ia /°K; A{T\ 



Since in both isothermal and thermogravimetric cases usually the same analytic process is 
involved in activation energy determination, i.e., the logarithm of some function of a rate is plotted 
against l/T(°K), any error resulting from ignoring temperature dependence in the determination 
of E should be about the same in both cases. If A is linearly dependent on temperature, the dif- 
ferentiation of the logarithm of eq (14) with respect to \jT results in 



d In k__ T _E 
R 



© 



(22) 



Thus, at 7"= 500 °K, the correction for linear temperature dependence of A on the experimentally 
obtained activation energy will be about 1 kcal/mol. 

In order to observe the effect of linear temperature dependence of the preexponential factor 
on the shape of a thermogram, Case IV, £ = 60,000 cal/mol, Ai//8= 1.374 X 10 13 (°K)" 2 , was cal- 
culated for various orders using pi(x) values obtained from eq (16). The integral and differential 
curves for first order kinetics may be compared with Case I, E — 60,000 cal/mol, A/P = 10 ie (°K) _1 
in figures 8 and 9. The difference is trivial for these conditons. The effect is in the direction of 
increased slope and is greater for higher orders. It will be significant only at very low activation 
energies. 

If a single Arrhenius expression is applied to situations involving simultaneous reactions, the 
parameters E and A well may exhibit temperature dependence. This aspect is discussed in the 
section on composite cases. 

3.4. Effect of Rate of Heating 

Analyses of the changes in thermogravimetric data brought about by variation of the rate of 
heating, /3, are the basis for the most powerful methods for the determination of kinetic parameters 
and these methods are discussed in subsequent sections. 

Figure 10 contains the integral curves for Case B with maximum isothermal rate at 50 percent 
completion for A/p 2 = 10 16 /°K, £ = 60,000 cal/mol and heating rates of 01 = 0.05, j8 2 = 0.10, and 
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FIGURE 10. Effect of heating rate on residual fraction 
versus temperature. 

Case B: /4 = 10 ,!S /9ec, £=60,000 cal/mol, £,=0.05, /8,=0.10, 3 =O.2O °K/*ec. 



FIGURE 11. Effect of heating rate on thermo gravimetric 

rate versus temperature. 
Case B: ^=10 15 /sec, £=60,000 cal/mol, /3,=0.05, 0,=O.1O, /3 3 =0.20 °K/sec. 



£j 3 = 0.20 °K/sec. Figure 11 contains the thermogravimetric rate curves for the same case and 
conditions. Not only are the curves shifted to higher temperatures by increased heating rates 
but they become less steep. 

3.5. Composite Cases 

The majority of experimental thermograms found in literature represent cases in which two 
or more volatilization reactions take place. If the values for the kinetic parameters, e.g., E and 
A, for each reaction are such that the regions of weight loss occur at separated temperature ranges 
then a stepwise thermogravimetric trace results. In these cases, each curve between successive 
horizontal portions may be treated separately. However, many cases are more complex. An 
example is a polymer which undergoes both depolymerization and side group splitting or ring 
formation, followed by depolymerization of the thermally more stable product polymer. 

Here we investigate only two simple cases: 

(1) Two Independent First Order Reactions 

If a fraction of the material, a, volatilizes by a first order reaction with Arrhenius parameters 
A\ and Eu and the remainder of the material volatilizes by an independent first order reaction with 
parameters A 2 and E 2 , the residual fraction will be given by 



\-C = a 



PR 



+ (!-«) 



A*E* 



( RTJ 



(23) 



(2) Two Competitive First Order Reactions 

If the total material may be volatilized by two alternative paths, each having a rate propor- 
tional to the first power of all remaining volatilizable material and if the respective Arrhenius 
parameters are A\ and E\ and Ai and Ei for the two paths, then the residual fraction will be 



\-C = e 






p {-Rt)- 



M 



-SI 

RTJ 



(24) 
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FIGURE 12. Effect of heating rate on. residual fraction 
versus temperature for composite cases. 
£, = 30,000 cal/mol; At =4.458 X 10 fl /sec 
£, = 60,000 cal/mol; /* 2 = 10 ,5 /sec 
/9= 0.001, 0.01, 0.1 and 1.0 °K/sec 

Ca9e (1), independent first order reactions 

Case (2), competitive first order reactions. 




FIGURE 13. Effect of rate of heating on the thermogravi- 
metric rate versus temperature for composite cases. 

Ei, E-i. A u and A t as in figure 12 

(a) Case (2) competitive first order reactions 

(b) Case (1) independent first order reactions. 



To observe the effect of variation of heating rate in these composite cases, integral and differ- 
ential curves were calculated for a = 0.5, £"1 = 30,000 cal/mol, E 2 = 60,000 cal/mol, ^ 2 =10 15 /sec 
and ^i = 4.458 X 10 6 /sec for heating rates of 1.0, 0.1, and 0.01, and 0.001 °K/sec. The value of 
A\ was selected so that for Case (1) at /3 = 0.1 °K/sec and C = 0.5, the first reaction contributes 
75 percent to 1-C and the second reaction, the remaining 25 percent. 

The integral curves are shown in figure 12 and the differential curves in figures 13a and 13b 
for the two cases. 

Such a thousand-fold variation in heating rate is not easily realizable experimentally, but was 
used in these theoretical calculations to illustrate as wide a variety of situations as possible. 

It is immediately apparent from inspection of figures 12 and 13 and eqs (23) and (24) that for 
these cases of two first order reactions with widely differing activation energies, the thermograms 
for Cases (1) and (2) exhibit quite different variations in structure upon changing the heating rate. 

Case (1), (independent reactions), at the highest heating rate (and temperature), (/3 = 1.0), 
gives the appearance of one simple thermogram as the two reactions broadly overlap one another. 
As the heating rate is lowered, (/3 = 0.1), the curve tends to flatten as the two reactions begin to 
separate. The decrease in the value of the thermogravimetric rate at the maximum with decreas- 
ing j8 is the opposite of simple reactions as seen in figure 11. 

At /3 = 0.01, the two reactions separate further and two peaks are obtained on the differential 
curve. It may be noted that at these conditions the "maximum" temperature for the higher 
activation energy reaction has been shifted five degrees by the perturbing influence of the low 
energy reaction. 

By )8 = 0.001, the low energy reaction, favored at the lower temperature, is almost completely 
separated from the high energy reaction. 

On the other hand, in Case (2) (competitive reactions), for the values of parameters considered, 
the low activation energy reaction takes over almost completely at low temperature and heating 
rates (/3 = 0.001 and /3 = 0.01). At the high heating rate, (j8 = 1.0), the low activation energy 
reaction dominates only during the first 30 percent of weight loss, with the high energy reaction 
considerably modifying the latter portion of the thermogram. Again, the inclusion of the high 
energy reaction causes a temporary increase in (dC/dT) max with increasing (3. In simple reactions 
the maximum rate decreases upon increasing /3 and r max , [see eq (21)]. 

Two competitive first order reactions with equal activation energies (£ , i=£ , 2 ) appear as a 
simple first order reaction no matter how great the difference in preexponential factors. This 
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is true of any set of competitive reactions of the same order. If they are of differing order, then 
they are, in theory, resolvable. 

The differentiation between two independent first order reactions with the same activation 
energies depends upon how greatly A\\A% deviates from unity. The resolution is relatively un- 
affected by changing the rate of heating. 

Thus, the technique of varying the rate of heating in composite cases not only often permits 
a separation of reactions but often gives information as to the nature of the competition between re- 
actions. This aspect is discussed further in succeeding sections on methods of kinetic analysis. 

The application of more complex empirical substitutes for the Arrhenius equation to experi- 
mental data involving composite reactions has been discussed by Farmer [16]. 

4. Critical Survey of Methods of Kinetic Analysis 

The methods of kinetic analysis of thermogravimetric data are divided into five categories 
to facilitate discussion and comparison. These are -(a) "Integral" methods utilizing weight loss 
versus temperature data directly, (b) "Differential" methods utilizing the rate of weight loss, 
(c) "Difference-Differential" methods involving differences in rate, (d) Methods specially appli- 
cable to initial rates, and (e) Nonlinear or cyclic heating rate methods. 



4.1. Integral Methods 

The fundamental criticism of integral methods as applied to a single thermogravimetric trace 
is that "best values" of A, E, and n are inevitably fitted to the data whether or not these parameters 
have any significance or even utility in the understanding of mechanism. Also to obtain the activa- 
tion energy the order must be known or vice versa. This latter problem has been met with varying 
success in several ways. 

As many of these methods involve an approximate integration of the exponential integral, 
they will be compared to eqs (13a, b, c, d) when applicable. 

The first serious theoretical treatment of thermogravimetric data was by van Krevelen et al., 
[17]. These authors approximated the exponential integral by making the substitution 



_E_ 
RT 



-(.-¥)*.(, 



0.368 7^' Rl - 



(25) 



where T max is the temperature at the maximum thermogravimetric rate. Thus, they obtained for 
eq (4) in logarithmic form 



In 



(l-C)'-"-I 
1-1 

(»»*1) 

lnln(l-C) 
(n=l)J 



= ln 



AR / 0.368 \ ifc _ 
fiE\T m J "-x + D 



+ 



^max+D 



InT 



(26) 



and tested for linearity of In T for various values of n. Reactions are insensitive to order at low 
conversion so in order to obtain n from this or similar methods linearity must be established at 
high conversion. However, linear plots over much of the conversion range were obtained from 
eq (26), tt=l, for Cases B and C with maximum isothermal rates. Thus, the method should be 
limited to cases in which a known isothermal order is followed. 
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Baur et al. [18], in 1955 applied a graphical solution of the exponential integral to the case of 
a linearly temperature dependent preexponential factor (eq (15)) for data on the high pressure 
oxidation of metals and obtained good agreement between experimental values and theoretical 
values based on a parabolic rate law. Smith and Arnoff [8] in 1958 calculated a table of p(x) 
values for x = 1 to 50 and used it to fit data on the thermodesorption of gases from solids to first 
and second order kinetics by a method of trial and error. 

Doyle [10], in 1961, presented the most formal method for the precise curve fitting of a single 
thermogram using calculated tables of log p(x). However, as in other such methods, its applica- 
bility is limited to cases of known isothermal order. As the method is quite tedious, Doyle devised 
an approximate method based on the first two terms of the asymptotic expansion (eq (13a)). 



Pa,2(x) = e-EIRT[-I f ) 



(27) 



For cases of unknown order, Doyle [10] suggests making calculations as early in the reaction 
as experimental accuracy allows. However, even using C=0.05, one obtains appreciable errors 
in the determined activation energy when this method is applied to maximum Cases A, B, and C. 
Nevertheless, such integral curve fitting of a single trace is more generally applicable than similar 
differential methods. 

Farmer [16] defines the relative error using the two term asymptotic approximation as 



r 2 z 



. P(x) 
Pa,2(x) 



(28) 



This relative error, r 2 , and r 3 , the relative error in using three terms of the asymptotic approximation, 
are shown for several values of x in table 1. 



TABLE 1. Relative errors for two- and three-term asymptotic ap- 
proximations o/p(x) a as a function o/x = (— E/RT) [75, 16] 


— x 




10 


20 


30 


40 


50 


r 3 


0.844 
1.055 


0.913 
1.014 


0.939 
1.006 


0.954 
1.004 


0.962 
1.002 



"PW = 



X J-oo X 



Farmer [16] substituted eqs (27) and (28) into eq (10) to obtain in logarithmic form 

f c dC_ 
Jo f(C) _ _E_ 



log : 



J»2 



1 

2.303/? \T 



-log 



or 



(29) 



C c dC 

Jo /(C). 



T log ^jp=-T log [r*j£[ + J. 



3037? 



where rz is a slowly changing function of temperature. The second expression was found to be 
sensitive to data point deviation. 

As r 3 lr 2 =l-\-2jx=l—(2RT/E), if the first three terms of the asymptotic approximation are 
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utilized and if /(C) is given by eq (3) then 

"1-(1-C)'-"1 



log 



TH 



log 



-C)'-» ] 
l-it) J 

(n*l) 

1-C 



In 



r 2 



slog 



^ 

££ 



1 2/m 



£ J 2.303/?r 



(30) 



which was developed by Coats and Redfern [2, 19]. A plot of log F(C)/T 2 versus 1/T appears to give 
quite satisfactory results when the kinetics follow a simple order. However, as with previous 
methods, application to systems of "changing order," e.g., maximum Case B, results in erroneous 
values for E. 

As Doyle has done, Coats and Redfern [20] have gone to low conversion data to get around the 
problem of determining order. 

Assuming all reactions behave as zero order as C— >0, they obtain from eq (30) 



. C , AR 
log ^ = log — 



2RT \ 



1 



2.303/? T 



(31) 



and a plot of log C/T 2 versus l/T for low C or n = should give a straight line with slope of — E/4t.6 
as log ARI(3E[l — 2RT/E] is again sensibly constant. It is doubtful that this approximation may 
be generally used up to 6" = 0.3 as is suggested as for at n = 2, the error in F(C) amounts to 40 
percent at C = 0.3. 

Horowitz and Metzger [21] simplify the exponential integral with an approximation similar 
to but simpler than van Krevelen et al., i.e., defining a characteristic temperature, 0, such that 
6=T—T S where T s is a reference temperature at which \ — C= lie. Then making the approxi- 
mation 






1 



1 



Ts + 



°K) 



T s ' 



e 



(32) 



they finally obtain for n= 1, 



lnln(l-C); 



E6 



(33) 



so that a plot of In ln(l — C) versus 6 should give linearity with a slope of E/RTs . 

For cases in which n ^ 1, more complicated expressions are suggested [21] involving derivative 
parameters and T s is defined as the temperature at the inflection point. 

One may obtain from eq (12) using eq (13d) for \ogp(x) and making the approximation of eq (32) 



- |n r . . AE . M 1.052 , 1.052E6 

In ln(l — C) = ln -r^ — 5.33 ^ 1 — PT7 2 

pti Is tils 



(34) 



which gives the same slope as eq (33). However, eq (34) does not include the condition that T s 
is the temperature at which 1 — C=l/e. 
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Thus, eq (34) suggests the generalization of Horowitz and Metzer's method to the case where 
n^l and Ts is any reference temperature as in eq (35). (The "constant," 1.052, may be improved 
upon once an approximate E has been determined from a table of first differences of log p{—EjRT) 
for various E/RT, [11].) 

ln[l-(l-C) 1 - w ] = constant + L ^ 2 f (9 (35) 

E may be calculated using eq (33) with fair agreement with theoretical values for E = 60,000 
cal/mol, A/(3= 10 16 /°K for ra=l, but theoretical data for n= 1/2 as well as, for example, maximum 
Case C, also give nearly linear plots for In ln(l — C) versus 6 over a wide range and high values 
for the energy of activation. 

Few papers on computational methods have been published [2, 22], but such methods are 
probably widely used. McCrackin [23], in an unpublished work, has programmed a method in 
which a series of weights, w\, and temperatures, 7*, are fed to the computer and values of (1 —C)i 
and F(Ci) are computed for each temperature, assuming n. From equations 



FiCi)= L z=cTJ* 



(36) 



where Xi = (EjR)T-\(EIRT) and T_i is the incomplete gamma function of— 1 order, values of xi are 
computed for an assumed E and data fitted to eq (36) to determine A//3. Assuming the error in 
d is only experimental and independent of its value, F(Ci) will have a constant variance so the 
best estimate for A/fi will be given by 

2xiF(Ci)(l-Ci) 2n 

i 

From eq (36), F(Ci) is computed for each i using the calculated values of A//3. Then, (1 — d) 
and W\ are computed. The standard error of fit is 

, ^ [r ? (calc)-rKexpt)] 2 
o-{n, E) = jy—^ (38) 

N = No. of values of T and W . 

This is computed for various values of n and E and the values that give the smallest cr are 
chosen. 

To summarize, the computer calculates and prints cr for values of E for n = 0, 1/2, 1, 3/2, 2, and 
3. Then, for each /i, it chooses the value of E that gives the minimum cr and prints E, cr and A\$ for 
each n. It then prints ^(calc) — Wt{expt) for these values of E. 

Since only if Wi{ca\e)— Jl^(expt) remains sensibly constant over the entire range can the system 
be said to have a particular n and E, it would appear that McCrackin's method tests the validity 
of its assumptions. 

Reich and Levi [24] have developed a method involving the determination of the area under 
the initial part of an integral thermogram to obtain an approximate expression for the energy of 
activation. Doyle [25] points out that such a method should be less sensitive to experimental 
scatter than line or slope methods. The corrected derivation [25] of the equations includes the 
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approximation in eq (13d), i.e., p(x) = Ge x , to obtain for n=l, 



. . . A {EGY E 



which for two different initial areas becomes 



In/4+u 







InC 






(ra=0) 


In 


In 


(1-0 



where A l + u IA l is the ratio of initial areas. As will be shown subsequently, the coefficient of E/RT 
assumed to be unity in eq (39) is a slowly changing function of E/RT and more accurate values may 
be obtained for E by successive approximations. Also, the general utility of eq (40) depends on 
the lack of sensitivity to order in the early conversion range. 

Reich [26] has pointed out that at small and constant AT from eq (11) 



s -§f +ln i AT (41) 



AT= constant. 

Utility of eq (41) is limited to cases in which a particular reaction order, /i, is known to be 
followed. 

Reich [27] using the approximations made by Horowitz and Metzger [21] obtained at constant 
weight loss for two different heating rates, 

4.6 log | (f 

/V\ (42) 

by assuming that T s , ilTs,2 = Ti/T 2 . 

However a simpler expression has been developed [11] by the substitution of eq (13d) into 
eq (10) to obtain 

AE E 

log F{C) = log — - log fi- 2.315-0.457 ^ (43 ) 

so at constant conversion for several heating rates a plot of log f3 versus 1JT will have a slope of 

^ - 0.457 | or E = - 4.35 ^ (1 - C) = constant. (4 4) 

d— d — 

Activation energies may be quite accurately and simply obtained from eq (44) by successive 
approximations as tables of log p{— E '/ RT) and A log p(—E/RT) for various E/R T are available [10, 11]. 
Once an approximate E is obtained from eq (44), the "constant," which changes from 0.477 at 
E/RT = 20 to 0.449 at E/RT= 60, is redetermined for the approximate E/RT and successively more 
accurate values of E are obtained. 
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The independence of E from C and T may be tested by determining E at various constant 
conversion values. 

In an excellent and important paper, recently to come to our attention, Ozawa [27a] also de- 
velops eq (44) but employs it without further refinement to calculate E at several C values. He 
sets up appropriate theoretical master-curves of (1 — C) versus log <$>[=log {AE/ftR)p{x)] for, e.g., 
simple nth order reactions and for the random degradation of polymers [eq (6)] for N>L and 
L = 2, 3, 4, and 5. A more accurate experimental master curve is obtained by superimposing 
curves of log ft versus l/T at several heating rates by displacement along the abscissa. Values 
of C and T from this curve are used to construct a plot of (1 — C) versus log [(E / 'ftR)p(E / RT)] which 
may be matched to the appropriate theoretical curve by displacement of log A along the abscissa, 
thus determining^ and confirming the assumed kinetic equation. 

This appears to be one of the best and most generally applicable methods yet developed. 

If F(C) may be represented as a simple nth order reaction as in eq (3) then eq (44) becomes, 



log 



(1-C) 1 -*-! 
l-n 



g ~R — og ^~ 



2.315-0.457 — 



(n*l) 



(45) 



so for the slopes of plots of log ft versus l/T at constant C, log [(1 — C) 1 n — 1] versus l/T at constant 
ft, and log [(1— C) 1_n — 1] versus log ft at constant T, one obtains 



d log ft 



E _d\n [(1-C) 1 -*-!] 



.1 - "•*"*- 


d Tjn 


(C = constant) 


{ft — constant) 


dlog [(1-C) 1 -*-!] 




d log ft 




in *1) 


= 1 


dlogln(l-C) 




dlog ft 




(n=l) 






(T= constant). 



(46) 



(n*l) 



(47) 



The right-hand side of eq (46) may be applied to a single thermogravimetric trace and has the 
advantage over the approximate methods of Farmer [16] and Coats and Redfern [19], [eq (29) and 
(30)], in that successive approximations in the manner described above can be used to determine 
accurate values of the parameters. The order, n, may be tested by eq (47), the correct n giving 
a horizontal plot. 

Once ft, n, E, T, and C are known, one may substitute into eq (45) to determine A. Thence, 
constancy of A with changing T and C should be a confirmation of the validity of n. 

As this method is the most adaptable of the integral methods, we demonstrate its utility with 
the theorietical cases of "changing order" (maximum Cases A, B, and C) and "changing activation 
energy" [composite Cases (1) and (2)]. 

Figure 10 contains the integral curves for /3i = 0.05, ft 2 = 0.10, and /3 3 = 0.20 °K/sec for Case 
B with a maximum isothermal rate at C = 0.5. Values of T at constant C were obtained from these 
curves for C = 0.02, 0.10, 0.25, 0.50, 0.75, 0.90, and 0.98 and plots of log ft versus l/T are shown 
in figure 14. The slopes of these lines give, from eq (44), ^approx. = 59,300 ± 200 cal/mol. Upon 
reevaluating the constant in eq (44) for (^//^T^approx.^ one obtains ^corrected = 59,850 cal/mol com- 
pared with the theoretical value of 60,000 [11]. 
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FIGURE 14. Logarithm of heating rate versus reciprocal 
absolute temperature for maximum Case B. 

A= l() ,5 /sec; £ = 60,000 cal/mol 
0i = O.O5,/3,=O.lO£s=»O.2O°K/sec; 

r; = 0.02, 0.10, 0.25, 0.50, 0.75, 0.90, and 0.98. 




1.40 1.60 

l/TXI0 3 (°K) H 



Figure 15. Logarithm of heating rate versus reciprocal 
absolute temperature for composite cases. 

E\, E>. A\. A>, and /3 as in figure 13 

(a) Case (2) competitive first order reactions 

(b) Case (1) independent first order reactions 

#. C = 0.10 □, C = 0.75 

V, C = 0.25 O, C = 0.90. 

A, C = 0.50 



Therefore, this not only permits determination of the correct activation energy for this case 
of "changing order" but establishes its independence of C and T over the reaction range. 

Equation (44) should not be expected to apply to composite cases since it was derived assuming 
a single temperature dependent rate constant. However, application to the independent and 
competitive first order reactions shown in figure 12 and 13 demonstrate some of the limitations of 
the method and how it may be utilized to interpret these simple cases. 

In figure 15, log fi is plotted against l/T for — (a) two first order competitive reaction, and 
(b) two first order independent reactions, for E x = 30,000, £2 = 60,000 cal/mol, fi = 0.001, 0.01, 0.1 
and 1.0 °K/sec and C = 0.10, 0.25, 0.50, 0.75, and 0.90. The theoretical slopes for E x = 30,000 and 
£2 — 60,000 cal/mol are included as guidelines. 

For two competitive reactions, (fig. 15a), at low conversion, (C = 0.10), the low activation energy 
reaction is dominant at all heating rates so a slope corresponding to 30,000 cal/mol is obtained. 
At high conversion, (C = 0.90), the low energy reaction is still dominant at slow heating rates but 
the intrusion of the high activation energy reaction can be observed at fast j8. A readjustment of 
the parameters could set up a situation in which the high energy reaction was dominant except 
at low C and /3. 

For two independent first order reactions (fig. 15b), at low conversion and slow /3, the low 
energy reaction is dominant but at low C and fast j8 the slope is perturbed by the high energy 
reaction. At high conversion, the high energy reaction is dominant at slow /3, but, as before, 
a mixture of £1 and E 2 reactions contribute at fast /3. 

In general, employing the lowest practicable heating rate will best isolate competing reactions. 
However, it may occasionally be expedient to raise the heating rate to pick up a high energy reac- 
tion that may not be easily discernible under near-isothermal conditions. 

We conclude, therefore, that only methods involving several heating rates can give the correct 
activation energy for cases of "varying order" (Cases A, B, and C) and reveal and, to some extent, 
resolve cases of "varying activation energy" (Cases (1) and (2)) in which several competing reactions 
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The integral methods involving a single thermogram appear to be applicable only in special 
cases in which the isothermal order is known and meaningful. However, integral curve fitting 
should be considered important corraborative evidence in the testing of mechanism. 

4.2. Differential Methods 

"Differential" methods based on rate of weight loss versus temperature data have been 
devised which are much simpler in application and, in some cases, are able to circumvent difficul- 
ties found in many "integral" methods. However, they suffer from an inherent weakness — the 
magnification of experimental scatter — often rendering their application to experimental data 
difficult, if not impossible. 

Van Krevelen et al. [17], calculated solutions for p(x) using two terms of eq (13c) and plotted 
families of curves for various T max for log {T(dCldT)) max and log (AT/T) max versus log E/R for 
first order reactions. AT is the half- width of the differential curve. 

Turner and Schnitzer [28] used three terms of eq (13c) for p(x) and refined van Krevelen's 
relationships for n = 1 to obtain expressions relating E/R to T max and A7\ A7\ !T max and (dCldT) max 
were estimated from a differential curve for each of several composite reactions and values for 
the parameters obtained. Turner et al. [29], using an identical development, obtained similar 
expressions for n = 2 /3. However, the equations in references [28] and [29] are in error [15]. 

Kaesche-Krischer and Heinrick [30] used van Krevelen's method at several heating rates to 
investigate the composite kinetics of polyvinyl alcohol decomposition. The utility of this method 
depends on how closely the assumed order of kinetics is followed. If the kinetics, in reality, 
behaves as a lower order reaction, then the calculated value for E will be too high, and conversely, 
if n is greater than its assumed value, the calculated E will be too low. Application of this method 
to maximum Cases A, B, and C gave activation energies from 70 to 150 kcal/mol compared with 
the theoretical value of 60. Van Krevelen [17], assuming n=l, obtained activation energies for 
polystyrene and polyethylene of 82 and 98 kcal/mol. These polymers which often exhibit maximum 
type isothermal rate curves usually have activation energies of the order of 55 and 70 kcal/mol, 
respectively [5]. 

Equation (11) at the maximum rate becomes 



//iC\ A 

=j (l -Os^e-*/" w (48) 



\dT) n 



which may be combined with eq (21) and eq (30) which utilizes the three term asymptotic approxi- 
mation to obtain 

±v± max H 



and 



"(l-C)-^! +(*-!) ^jp (50) 



In (1-C) max = 



(rc>0, 7^1) 

2/?r max i 

E J 
(n=l). 



These eqs (48), (49), and (50) were first derived for first order reactions and used by Murry 
and White [31]. The solution of eq (49) was facilitated by using several heating rates. Kissinger 
[32] extended Murry and White's method to any order, n, by deriving eqs (48), (49), and (50) as 
shown above. 
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Kissinger [33] differentiated the logarithm of eq (49) to obtain for n = \ 



din 



T 2 F 

1 max ^_ *-* 



R 



(51) 



Thus the activation energy may be obtained from the shift of the maximum temperature with 
heating rate. Equation (51) is identical to eq (42) which was also derived assuming n=l and a 
reference temperature [27, 21]. Substituting eq (50) into (49) and differentiating the logarithm 
assuming T max changes slowly with /3, Kissinger [32] again obtains eq (51) now for any n. 

Kissinger [32] also has developed a shape index, s, defined as the absolute value of the ratio 
of the tangents to the differential curves at the inflection points and related to the order by 



n=l.26s 1 l 2 . 



(52) 



Such higher derivative methods are seldom practical for polymer degradation studies. 

Fuoss et al. [34], suggest determining the three maximum values, T max , (dCldT) max and (1— C) max 
from the inflection point of the integral curve and determining the activation energy from eq (21). 
The method as described [34] is applicable only to first order kinetics and not to all orders as is 
implied by the omission of n from eq (21). Such an omission leads, obviously, to the exceptionally 
high value for E for polystyrene indicated in this paper [34]. 

However, (1 — C) max is relatively independent of the heating rate, /3. If eq (10) including two 
terms of the asymptotic approximation, eq (13a), is combined with eq (49) at the maximum, one 
obtains 



(l-Omax^rt 1 



(/l>0, 7^1). 



(53) 



Equation (53) was first pointed out by Kissinger [32] and also was derived by Horowitz and 
Metzger [21] from their nearly equivalent approximation. Thus eq (21) becomes 



11 IX1 max 



(tt>0, #1). 



(54) 



The last two columns in table 2 give the approximate (1 — C) max , [i.e., columns (1 — C) max {x = °°) 
and n n l n ~ 1 R(x = <*>)], which may be utilized with eq (54) to determine activation energies if the 
value and constancy of n have been established. If (1 — C) max is independent of /3, then from 
eq (21) T^^idC I dT) max must also be independent. Such independence is inherent in the use of 
the two term asymptotic approximation as differentiation of eq (48) with respect to /3 assuming 
T^^dCjdT) max and (1 — C) max constant results in eq (51). 



TABLE 2. Values of {1 — C) ma x for various orders and x's 
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1/2 


0.294 


0.273 


0.266 


0.262 


0.260 


0.250 


3.97 


1 


.430 


.401 


.391 


.385 


.382 


.368, l/e 


5.40 


3/2 


.516 


.484 


.472 


.465 


.461 


.444 


6.71 


2 


.578 


.544 


.531 


.523 


.519 


.500 


7.95 


3 


.661 


.626 


.612 


.603 


.599 


.577 


10.32 


4 


.716 


.681 


.666 


.658 


.653 


.630 


12.62 


5 


.755 


.721 


.706 


.698 


.693 


.669 


14.85 
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It also follows that 



cy CI ill 1 max __ \Cll } max ___ 1 (ZEi\ 
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and from eq (8) that 
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so, if n is known, £ and A may be determined from a plot of log /-. _ /^\™ aX versus ^ — for several 

I* ^Anax * max 

/3. Differentiation of eq (56) with respect to j8 assuming (1 — C) max constant results in 
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(I), 



R (57) 



which is an equivalent form of eq (51). 

Farmer [16] has developed equations for 7*1/2, the temperature at which C = 0.5. The relative 
error [eq (28)], r2, in using the first two terms of the asymptotic approximation is related to con- 
version at the maximum by [16] 

(l-C) max =e- r2 

(»=1) 

i 



(l-Omax=[4 ZL r 2 + l 



(58) 



Therefore the two-term approximation holds exactly (r 2 = l) only at x=—EjRT=^. Utilizing 
values of n> as functions of x from table 1, we can calculate the (1 — C) max values for several values 
of x as is shown in table 2. 

The approximate values for (1 — C) max (— # = «>) may be in considerable error for cases of low 
EIRT max and use of eq (53) will result in large error in determining n. However, if the approximate 
E/RT max is known, appropriate corrections may be made from table 2. 

If dC/dl/T is plotted versus l/7\ the maximum rate is 



•(§. 



\E *" ^^max) yl ^)max 

nR (21a) 



which is practically independent of temperature. Other equations at the maximum will be modi- 
fied accordingly. 
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The actual application of these many mathematical excercises at the maximum rate are limited 
to cases where n is known and its validity over the entire range of conversion has been established. 
An example of the magnitude of the error which may result in the misapplication of a single point 
method may be observed from the calculation of E from Case B. First order kinetics would appear 
to hold as (1 — Q max =0.35 for this case. However, calculation of E from eq (21) for n = \ gives 
^apparent — 135,000 cal/mol compared with ^theoretical = 60,000 cal/mol. Subsequently, it will be 
shown that a first order curve will give a deceptively good fit for both the differential and integral 
curves for such cases of maximum isothermal rate. 

The simplest differential method for determining kinetic parameters is an Arrhenius plot of 
the logarithm of the rath order rate constant, k n , against the reciprocal of the absolute temperature 
since from eq (11), 



log k n =log 



or at j8 = constant, 



P 



dC 
dT 



(1-Q» 



= logj8 



rfln(l-C) 
dT 



-{n-\)\og{l-C)=\ogA- 



2.3RT 



(59) 



log 



dC -, 
dT 

d-cr 



= log 



fi 2.3RT 



Kofstand [35] used eq (59) to test for linear (n = 0), parabolic (n=— 1) and cubic (n=—2) rate laws 
and determine activation energies for the oxidation of metals. 

Barrer and Bratt, Newkirk, and others [36-40, 2] have employed eq (59) for n=\ to determine 
activation energies from thermogravimetric data. This method may be used early in the reaction 
to determine initial parameters as the early parts of most reactions are more or less independent of 
order. However, it has been suggested that this method be used over the whole range of volatiliza- 
tion [2,37] with values for E obtained from linear portions of the log (din C)/dT versus 1/Tplot. As 
with other differential procedures this will lead to too high values for E if n < 1 and too low values 
if n> 1. Thus, when this method is applied to Cases B and C, high values are obtained for E. 
In order to obtain accurate activation energies, one would have to know or guess the approximate 
value of n that the reaction appeared to follow in the range of determination. 

Ingraham and Marier [41] have used this method for a zero order reaction assuming a linear 
temperature dependence of the preexponential factor and obtained a linear relationship between 
log UT(dCldT) and l/T. 

The differential methods treated thus far assume the existence of a single order, n. At worst, 
they determine order at a single point such as the maximum rate. At best, they do not rigorously 
and sufficiently test their initial postulate as to the constancy of the parameters as does the inte- 
gral method at several heating rates or as does a similar differential method which follows. 

As in the case of integral methods, it is necessary to perform experiments at several heating 
rates in order to reveal changes in the kinetic mechanism which might affect the reliability of param- 
eters determined from the data. Anderson [42] solves the three simultaneous equations for eq 
(59) at three different ft with an electronic computer for the parameters A, n, and E at a series of 
constant (1 — C) values. Any systematic drift in parameters with changing conversion would be 
indicative that the assumptions — single order and single temperature dependent rate constant — 
do not hold. 

The most generally applicable differential method was developed by Friedman [43] who utilized 
eq (9) in logarithmic form, i.e., 



In 



%-*fig-*4fiC>-£r 



(60) 
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{dC/dt) and T are determined at constant C for several /3, thus a plot of In (dC/dt) versus l/T will 
have a slope of EjR and an intercept of In Af{C). If E proves to be constant over a range of C 
values then, from eq (59) and (60), a plot of In Af(C) versus In (1 — C) will have a slope of n and 
intercept of In A if an nth order kinetic mechanism is operative. 

In general, differentiation of eq (60) at constant C or constant dCjdt gives 



,, dC 

dt _E_ d\nAf(C) 



i) R 



d\ 



(61) 



orif/(C) = (l-0», 



(!) 

in * ( dC \ 

(C =const) I -j— =const I 



am dt _ E_ din (1-0 
1 — ~R~ n i 

^ d^ (62) 

(C = const) I —r-= const I 



, 1 dC 

a In -7- 

dt 

d\n(l-C)~ 
{T= const). 



(63) 



The validity of /(C) = (1 — C)* 1 may be tested from several runs at different (3 either at constant 
rates or constant temperature. 

The right-hand side of eq (62) may be applied to a single differential plot when n > 0, since 
under these conditions dCjdt is a dual-valued function. The constancy of E/n may be tested at 
a series of constant dC/dt values. However, such constancy does not "prove" that the reaction 
is nth order. Experiments at several heating rates are essential for such proof. 

Chatterjee [44] quite recently has expressed eq (63) in terms of W, the weight of sample 
remaining, 

(T= const) 

and the order n may be determined at several temperatures from two runs differing only in their 
initial weights. 

It may be noted that the preexponential factor of the dC/dt equation is JF8" 1 times the corre- 
sponding preexponential factor of the dW\dt equation [45]. Also, if an initial order is determined by 
this method, a distinction between the "order with respect to initial weight" and the "order with 
respect to temperature (time)" may be useful in the interpretation of the mechanism of complex 
reactions. The order with respect to initial weight, n w , may be determined from two or more runs 
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at different Wo from an integral method, since from eq (43) 

Alog)To _ 0.457£ 

4 ' 



(n w -l)R 



at constant W and /3. 

After the determination of n, Chatterjee [44] suggests determining E from a single run by a 
difference -differential method. However, his method may be extended to apply to any j\W) ¥* W n 
since eq (62), for two or more runs at different initial sample weights becomes 



din 



dW 

dt _E_ dlriAf(W ) 



,1 

Am 



R 



1 
T 

dft 



(65) 



(const W) I const -t— J 



Chatterjee [44] suggests using eq (64) to determine the order of competitive reactions by first 
calculating A\ and E\ early in the first reaction and subtracting off the portion of the rate (calculated 
from A\ and E x ) due to the first reaction from the second and so on. When change in n is due to 
competitive reactions this method may be applied only when the orders of the several reactions 
are different [16] and the rates of each step are proportional to a power of the total weight of mate- 
rial. An added problem is the determination of n early in the reaction. 

We may test the applicability of Friedman's method [43] to the cases of "changing order" 
and "changing activation energy." Figure 16 is a plot of log dC/dt versus 1/7" for Case B with 
a maximum isothermal rate at 50 percent conversion. As with the integral method, parallel straight 
lines are obtained for degrees of conversion ranging between 0.04 and 0.95. The average value 
of £ ca . c is 59,800 ± 1500 compared with 60,000 cal/mol, theoretical. 

In /3 (dC/dT) is plotted against 1/7 1 lor two first order independent composite reactions in 
figure 17, as was done for the integral method in figure 15b. Figure 18 is a similar plot for two 
first order competitive reactions. Again, theoretical slopes for E x =30,000 and E> =60,000 cal/mol 
are included for comparison purposes. The close similarity of these results to those of the integral 
method is obvious and the interpretation of figures 17 and 18 is identical to that of figure 15. Com- 
parison of eqs (61) and (44) shows the dominance of A In f3 over A In [dCjdT) and thus, explains 
the close similarity of the two methods. 



Figure 16. Logarithm of thermogravimetric rate versus 
reciprocal absolute temperature for Case B. 

/. = 10 I! 7sec; £ = 60,000 cal/mol 

/3, =0.05. j8-. = 0.10. #, = 0.20 °K/sec 

C = 0.04, 0.10, 0.20, 0.30, 0.40, 0.50, 0.60, 0.70, 0.80. 0.90, and 0.95. 
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Figure 17. Logarithm of thermo gravimetric rate versus 
reciprocal absolute temperature for Case (1) independent 
first order reactions. 



FIGURE 18. Logarithm of thermo gravimetric rate versus 
reciprocal absolute temperature for Case (2) competive 
first order reactions. 

£,, E>. A,. A-,. (3. as in figure 13 



•. C = 0.10 


• D, C = 0.75 


V. C=0.25 


O. C = 0.90 


A. C = 0.50 





E^. 


E, 


A t ,A t 


/3 


as in figure 13 


• 


€= 


= 0.10 




□, 


C = 0.75 


V, 


C- 


= 0.25 




o. 


C = 0.90. 


A 


C- 


= 0.50 









The practical merits of the two methods may now be compared. The integral method [eq 
(44)] is simpler than the differential [eq (61)] method and quicker as it does not involve the deter- 
mination of rates. However, no successive approximations are necessary for the differential 
method. (It will be shown that the integral method may become exact by appropriate programming 
of the heating rate.) 

On the other hand, many polymer decompositions are subject to early kinetic irregulari- 
ties (e.g., a temperature dependent induction period). Such complications may modify param- 
eters obtained from an integral method as they depend upon cumulative values of the experimental 
parameters. Differential methods which give instanteous values for parameters are not subject 
to these complications. The preexponential factor, A, and the order, n, if it exists, may be deter- 
mined directly from the intercept in either method and both methods are equally capable of inter- 
preting cases of changing order or activation energy. 

To summarize, of the differential methods, only those involving several thermograms appear 
to be generally applicable. Methods involving a single thermogravimetric trace should be applied 
only to systems where all material volatilizes by the same simple kinetic process. 

4.3. Difference-Differential Methods 

The difference-differential method of Freeman and Carrol [46] and its modification [38] is 
the most widely used method for the kinetic analysis of thermogravimetric data. It has been 
applied both to the investigation of inorganic materials [46, 2, 39, 47, 48] and polymers [38, 49, 2, 
50, 44]. 

The difference form of eq (59) is 



dC . _ „dC . _ „ ^ E 4 1 
Alog — Alog/3-^AlogXl-O-^A- 



(66) 



and the difference forms of eqs (62) and (63) are 

dC 



Alog/3 



dT 



Alog(l-C) 



»(?. 



2.3/?Alog(l-C) 
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(67) 



Figure 19. 



Mog 



dC 
dt 



Mog(l-C)] 

/4=10 15 see; E 

intercept — n 



lMog(\-C)j 

60,000 cal/mol; j8= 0.10 °K/sec 
slope =E/2.3R. 



for Case B. 




2.3/i# A log (1-C) 2.3# 



A 1 r> ^ 

A log £ ^ 



AlogiS 






(68) 



A log j8 



^/r_ Aiogd-o 



© *G 



2.3/2 



(69) 



Equations (67) and (69) have been used to obtain the parameters, E and n, from thermogravi- 
metric data with reported success not only for simple inorganic decompositions but also for com- 
plex polymeric systems. Therefore it is of interest to apply these equations to Case B with a 

dC 



maximum isothermal rate at 50 percent conversion. Figure 19 is a plot of 

x 1 



A log (1-C) 



against 



— — — for this case with a constant energy of activation of 60,000 cal/mol. Similar results 

A log (1-6) 

are obtained from eqs (68) and (69). 

Three conversion ranges of constant slope may be obtained from figure 19 as might be antici- 
pated for a cubic. These give the following values for n and E, 



1-13% conversion 
13-50% conversion 
50-95% conversion 



E= 66,500 cal/mol 
£=105,000 cal/mol 
E= 175,000 cal/mol 



n = -2A 
rc = + 0.3 
n = 1.25. 



The data may be interpeted as follows: between 1 and 13 percent conversion, the data are 
fitted by the parameters £ = 66,500 cal/mol, and n = — 2.4, etc. Since n is negative initially and 
becomes positive the isothermal rate must pass through a maximum value. The value of E does 
not have any correspondence to its theoretical value except in the initial conversion range. Treat- 
ing each linear range independently does not improve the results. Therefore, the difference- 
differential method gives a procedural n and E for the different ranges but these have no mech- 
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Figure 20. 



Mog 



dC 
dt 



A M1/T)J 
petitive first order reactions. 
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Figure 21. 
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Eu E->, At, A->, and /3 as in figure 13. 



for com- 



anistic significance, e.g., the activation energies obtained for the middle conversion range of 
polystyrene and polyethylene in reference [38] may well be in error. 

The Freeman-Carrol method is applied to two competitive first order reactions (Case 2) of 
30,000 and 60,000 cal/mol, respectively, in figures 20 and 21. 

Analysis of these figures allows a critical investigation of this method as first order kinetics 
is followed by both reactions. The scatter of points at low conversion in figure 19 due to the 
determination of 1 — C and dC/dt values from theoretical curves was cut down by determining them 
directly from the theoretical equations for figures 20 and 21. 



A log 



d£ 

dt . , , . A log (1 
- is plotted against — - 

A if 



-Q 



in figure 20. The theoretical intercepts at C— » for £" = 30,000 and 60,000 cal/mol and the the- 
oretical slopes for n—\ are indicated. As C-^>0, the curve approaches the correct intercept, 
(E= 30,000), but the slight perturbation by the £ = 60,000 reaction coupled with the lack of sensi- 
tivity toward n at low conversion tends to make the reaction appear to be zero order. The slope 
approaches n = \ only at high conversion, but with an intercept of £ = 40,000. 



A log 



dC 
dt 



A log (l-C) 



is plotted against 



4 



Alog(l-C) 



in figure 21. The theoretical intercepts at C— > 1.0 for n= 1 and the theoretical slopes for £=30,000 
and 60,000 cal/mol are indicated. In the less than 30 percent conversion range, a slope reasonably 
close to 30,000 may be obtained but the perturbation of the high energy reaction is great enough 
to throw the intercept toward zero order. The intermediate conversion range gives slopes between 
30,000 and 60,000 and intercepts less than unity. Only as C— > 1 may n be accurately determined 
from the intercept. 
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It is obvious from the above cases that it is difficult to differentiate between deviations from 
constancy due to "changing order" and "changing activation energy" with this method. 

If the reaction follows a simple order, // may best be determined at high conversion either from 
the slope of eq (69) or the intercept of eq (67). The activation energy may be determined either 
from the intercept at low conversion for eq (69) or the slope of eq (67). 

In general one may obtain the initial parameters only if values of high accuracy can be obtained 
for A dC/dt at low conversions. The experimental scatter, doubly magnified by taking the differ- 
ence of a derivative, often will not allow the determination of order at low conversion where the 
dependence on n is slight. 

These methods seem to be of limited applicability to polymeric systems where the kinetics 
often is complex. 



4.4. Initial Rates 

Even where /(C) in eq (1) is unknown, the activation energy often can be determined as near 
to initial conditions as possible since all well-behaved reactions approach zero order as C — » 0. 

Since the initial conditions for a thermally decomposing system are usually those most pre- 
cisely characterized, the initial values of the kinetic parameters are often particularly useful 
in determining mechanism. Accurate initial weight-loss values are unattainable in isothermal 
studies due to an inevitable time-lag in reaching experimental temperature. 

This independence of order at low conversion means that a plot of the logarithm of a "zero 
order rate constant," e.g., dC/dt, C/T 2 [20] , C [26] , or the logarithm of a "first order rate constant," 
e.g., (d In C)/dT [36-40, 4, 16], log (1 — C) [26], or the logarithm of an "nth order rate constant," 
[1/(1 — C)"] dC/dT, [35], versus 1/7 1 or l/T+a In T will give approximately the same values for 
activation energy at low C. Difference-differential equations, e.g., eq (69), to determine initial 
E have been suggested [44] but are not usually practical for reasons stressed in the previous section. 

The In T term which modifies the reciprocal temperature may enter as a temperature depend- 
ence of the preexponential factor, through an integral series approximation or a temperature 
dependence of the heating rate. It has been shown to be trivial except at low E. In general, 
if the order assumed in calculating E is higher than is actually the case, the resultant activation 
energy will be greater than the true value, although, where integral expressions are used other 
approximations may reverse this effect. 

From eqs (19) and (20) for n = 0, the slope of a plot of dC/dT versus C is EjRT 1 . [In fact, for 
any n, plots of [1/(1 - C) n ] dC/dT versus (1 - Cy- n ll~n, (n ^ 1), and d In (1-Q/dT versus In (1~C), 
(n= 1), have a slope of E/RT 2 .] E may be determined from the initial slope of a dC/dT versus C 
curve as the second term in eq (19) isnegligible at low C. This is shown for Case B in figure 22. 
From the slope, one obtains E^^^RT 2 (slope) = 670 2 X 2 X 0.0068 = 61,000 cal/mol compared with 
^theor — 60,000. The error in the calculated E results from difficulty in assigning a value to T, the 
average temperature over which the slope was determined. 

This problem may be overcome by plotting — dC/d(l/T) (= T^idCjdT)) versus C. The slope then 
becomes 



d / T2 dC\_ d 2 C 



dC \ dTj ir , . 1 

dtd — 



E 


+ 2T- 


f'(C) dC 


R 


AQ ,1 


E 
R 


+ 2T+ " dC 
1-C , 1 



(70) 
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FIGURE 22. Initial thermo gravimetric rate versus con- 
version for Case B. 

/*/j8=10"»/°K; £, he „ r = 60,000 cal/mol. 



Thus at low C, the slope will equal E/R + 2f where £//? ^> 2T. In practice an average E may be 
obtained from a set of C, dCjdT and 7 1 at low C since 



A [P 



dTj E — \ \dTJj _ r ^ 
A^ = R +2r ° r — C~- = E+2T ' 



RTf (§). _ 
£s ^LLi-ZRTi at C« 1. 



(71) 



Several methods have been applied to determine the initial activation energies for simple 
orders and the "maximum" cases. At C = 0.05, Doyle's method [10] give £ = 62,400 cal/mol for 
Case B. Equation (31) of Coats and Redfern [20] gave £" = 61,800 at C=0.04. Equation (41) of 
Reich [26] gave a low value of 59,000 as eq (41) differs from eq (31) by 2 In T. 

In general, (dCjdT) versus C plots remain linear to higher conversions than log (dC/dT) versus 
1/7 1 plots. Methods involving log k n versus 1/7 1 gave the correct slope up to C = 0.01 for Case B, 
while dC/dT versus C is linear to C = 0.03. 

If a procedural n may be approximated, log [(1/(1 — C) n ]dC/dT versus 1/r or [(1/(1 — C) n ]dC/dT 
versus (1 — C) 1_w /1 — n may remain linear over a wide range of concentration. However, a linear 
plot of log k n versus 1/7 1 at low conversion does not imply nth order kinetics as has been stated by 
some investigators. 

The apparent advantage of thermogravimetric over isothermal methods for the determina- 
tion of initial parameters is not always experimentally realizable. Accurate data are necessary 
at low conversions and should be obtained at low heating rates to separate any competing reac- 
tions such as volatilization of solvent., 



4.5. Nonlinear and Cyclic Heating Rates 

Techniques involving the programming of reaction variables are as yet a quite undeveloped 
tool for the elucidation of kinetic mechanisms. In general, we may write, 



dC = ldYj_ 
dt y\ dt 



~-f(YuY 2 , 



.,Y U . . .,rj 



(72) 
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where the reaction variables, F,, may include concentration of any reacting species, temperature, 
light intensity, dose of ionizing radiation, catalyst concentration, inhibitor concentration, product 
concentration, intermediate concentration, solvent concentration, dielectric constant of solvent, 
pH, ionic strength, etc. 

Time, of course, is a unique monotonously increasing variable. However all other reaction 
variables may, at least in theory, be programmed to increase or decrease with time during the 
course of the reaction at a rate given by 



dYj_ 
dt 



! A = *(2 1 ' i )i < 73 > 



f3 may be any function from zero (isothermal case, if Yj = T) to an integral transform and may in- 
clude several arbitrary constants and reaction variables. 

It may be simpler in some cases, e.g., the pH of an unbuffered solution, to allow a reaction 
variable to change in a prescribed manner rather than to keep it constant. If the change in a 
reaction variable may be measured, then an empirical equation may be fitted to its change with 
time. In such a case, any dependence of the empirical constants on the reaction variables should 
be investigated. 

Combining eqs (72) and (73), 

dc idc yif & Y ) 

dYi p, dt 



<2") 



A judicious selection of g may simplify the kinetics and information concerning the kinetic mech- 
anism may be elicited through the variation of the parameters. Also, this permits the analyses 
of a large variety of differential expressions of the form — dYj/dYj. 

An application of such techniques has been made recently by Ozawa [27a] and is well adapted 
to his method of kinetic analysis. For simple random degradation of polymers, the rate of change 
of fraction of bonds broken, da/dt, is given by a simple first order equation 

n da da . „,„„,. 

**'*-* {1 - al (75) 

Equation (75) may be combined, upon integration, with eq (5) to construct master curves of (1 — C) 
versus 4> for N>L and t = 2, 3, 4, 5, .... In general, if W=Wd\-C) = g{Y i ), and dYildt 
= Ae~ ElRr h(Y\), then integration results in relationship between conversion and O. 

Since E and the reduced time, d = <&/A= I e' E/RT dt= 0(/3, E, R, 7), are independent of C, 

J t 

Ozawa [27a] suggests that they be used to define thermal stability for materials for which the 
temperature dependence of the measured property can be expressed by a single Arrhenius 
equation. 

Many of the analytical kinetic procedures discussed in this article may be greatly simplified 
by the appropriate programming of the heating rate, for if 

^ dt 

and the temperature dependence of the preexponential factor is T\ T 1/2 , or T\ respectively and 
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then if m = 2, 2V2, or 3, respectively, one obtains 

1 dC dC A 



f{C)T~*e- 



P dt dT 
which is easily integrated to obtain 

J« f(C) aE J Xo a 



EIHT 
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-EjRT 



(76) 



(77) 



eliminating all problems associated with the exponential integral. 

Programming of /3 = aT 2 , i.e., time proportional to 1/T, should pose no special problem; for 
example, between 200 and 600 °C only a three-fold increase in heating rate would take place. 

All integral approximations are now unnecessary as, 
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and, for example, 
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The slope of the rate versus conversion curve now becomes, 

d 2 C = _E 2 /'(C) dC 
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The relationships at the maximum rate are simplified in some cases, 

. (E — 2RT max )(\ — C) max 
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or 

(dC\ E 
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nR°- C) - (21a '> 
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n[l-qti = l + (»-l)^|m ( 5o-) 



d In a £ 



\ i / max 



(51') 



etc. 

All differential and difference-differential equations are the same if /3(dCldT) = dCldt is 
replaced by aT 2 (dC/dT), thus, 
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ama dTE dln(l-C) 
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(61') 



etc. (62' 



and eq (63) does not change as is true of all equations at constant temperature. 

Thus, in summary, the programming of the heating rate proportional to the square of the 
temperature greatly simplifies integral methods rendering approximate equations exact. Rate 
versus composition equations are simpler as are expressions relating rate, composition and tem- 
perature at the maximum rate. Rate versus temperature methods are modified by the substitution 
dCldt = /3(dCldT) = aT 2 (dCldT) = -adCI(d IjT). 

Other cases of nonlinear heating rates have been developed for the investigation of the ageing 
of insulating material [51,52]. The case in which the rate of heating increases and decreases as 
an exponential function of time is of some interest as such behavior approximates the heating and 
cooling curves of an electrical oven heated at a constant load and then shut off to cool. Thus, at 
least in theory, once the time constant has been determined, such an apparatus may be used in 
thermogravimetric investigations without temperature control. 

The differential eq (9) is independent of heating rate so in many differential methods in which 
the simultaneous determination of T, 1— C, and dC/dt (dW/dt) are utilized, the analytic form or 
the value of the heating rate need not be known. Thus in eq (59), since /3 (dC/dT) = dC/dt, there 
is no need to know fi or to keep it constant unless a method comparing several heating rates is 
to be used. 

On the other hand, all integrated expressions assume a particular function for /3 so methods 
including integral forms depend on assumptions made concerning the rate of heating. 

Equations at the maximum, e.g., eq (21), are independent of the conditions used in reaching 
the maximum but assume that ft is constant in the region of the maximum, otherwise an additional 
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term containing the temperature acceleration would be necessary. Many of the other equations 
at the maximum are combinations of differential and integrated equations. 

Equations (61), (62), and (63) should be applicable in a two-point difference form irrespective 
of the previous variation of j3 for at any two temperatures at which (dC/dt) 1 = (dC/dt)2. 

| (jr- ~) = n In [(1 - C),/(l - C) 2 ] (62") 

or if the same temperature is obtained by any rate of heating followed by any rate of cooling, then 
at Ti = T 2 , 



ln [(f)/(f) 2 



ln[(l - C)i/(1 -C) 2 ]=n. (63') 



In order to utilize the left-hand sides of eqs (61) and (62), in a similar manner, two separate 
runs must be performed at different rates of heating such that at the same degree of conversion 
there will be different rates and temperatures so that at (1— C)i = (l — C)2 

Kf),/(f)L/GHH 

for any /(C). 

Reich, Lee and Levi [53] have applied eq (69) at constant dC/dT, 

A log _ A log W E 

~l~- n_ 7i Ur (69) 

A^7 A^; 

(dC \ 

\-^= const ) 

to cases where )8 was changed in magnitude and/or sign during a single run to obtain values for 
n and E. In runs in which the sample was both heated and cooled, eq (63') was also used to 
determine order. 

The above methods should be used with caution for polymeric materials. Except for the 
most simple systems, the initial starting material may undergo one or more condensed phase 
physical or chemical transformations before volatilization. In such cases, where two samples at 
the same degree of conversion have different heating histories, the rates of volatilization may not 
be comparable. Anderson [42] suggests, as a test for compositional constancy, that a series of 
M identical heating cycles be performed. If the relationship in eq (78) 

(1-C) M = (1-C)f (78) 

holds, then the kinetics follow a simple nth order relationship and a plot of the logarithm of the 
residual fraction after M cycles versus the number of cycles, Af, should give a straight line whose 
slope should equal the logarithm of the residual fraction after the first cycle. 

Anderson [54] also suggests that since for a case following a simple order n, 

f) M = M( l- W (79) 
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a plot of (dCldt)si versus M could be matched with curves for various orders. 
of \og(dCjdt)M versus M log(l — C)\ should yield a slope equal to n. 



Better still, a plot 



5. Concerning the Uniqueness of TGA Plots With Respect to Kinetic Parameters 

It is quite clear from previous sections that, even for cases in which a single Arrhenius expres- 
sion is operable, only methods involving several rates of heating determine a reliable activation 
energy parameter. Single point or even single curve methods may give spurious values for the 
parameters E and n as both the slope and the inflection point are dependent on /(C) and EjRT. 
A number of such methods give a reasonable fit for Case B (£ theor = 60,000 cal/mol, isothermal 
maximum rate at C=0.5) with a first order reaction but with an activation energy of approximately 
135,000 cal/mol. 

Therefore, in the final three figures, we compare these two cases in which the isothermal 
kinetics differ so drastically and observe to what extent their thermogravimetric curves are 
distinguishable. 

Three integral curves (1 — C versus T) are shown in figure 23. Curve 1 is for n— 1, £" = 80,000 
cal/mol and ,4//3 = 9.3179 X 10 21 /sec; Curve 2, Case B (isothermal rate rises to a maximum three 
times the initial rate at 50 percent conversion), £ = 40,000 cal/mol, Alfi = 1.2057 X 10 10 /sec, and 
Curve 3, Case B, £ = 35,000 cal/mol, A/fi = 3.6000 X 10 8 /sec. A//3 values were adjusted so that 
the temperatures of the maximum rates would coincide. 

The "maximum" Curves 2 and 3 deviate from the first order curves only in the first 20 percent 
conversion range. Better fit in this range could undoubtedly be obtained with a different value 
of n. 

Thermogravimetric rate versus temperature curves are shown in figure 24 for the same cases. 
The closeness of fit well explains the especially poor results obtained when methods involving 
maximum temperature, rate, conversion, or half-width are applied to Case B. 

Thermogravimetric rate versus conversion curves are shown in figure 25 for the same cases. 
The superiority of this type plot for distinguishing between the two cases is evident. The depend- 
ence of the initial slope upon only E and T permits a quick estimate of the initial activation energy 
whose constancy with increasing conversion may be established through other methods. 



■"========g g^rr---^^ 2 

V^3 




FIGURE 23. Residual fraction versus temperature. 



2. Case B: £=40,000 cal/mol; |= 1 .2057 X 10 w /«ee. 



FIGURE 24. Thermogravimetric rate versus temperature. 

I., 2., and 3. as in figure 23. 
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FIGURE 25. Thermogravimetric rate versus 
conversion. 

1., 2., and 3. as in figure 23. 



gree of 




To conclude, although there are some special cases for which a theoretical order has real 
significance, in general, n must be looked upon as a purely empirical parameter, sometimes useful 
in curve fitting. 



6. References 



[1] J. H. Flynn and L. A. Wall, Am. Chem. Soc. Meeting, Polymer preprints 6, 945 (1965). 

[2] A. W. Coats and P. J. Redfern, Analyst 88, 906 (1963). 

[3] C. B. Murphy, Analytical Chemistry 36, 347R (1964). 

[4] D. A. Smith, Rubber Chemistry and Technology 37, 937 (1964). 

[5] L. A. Wall and J. H. Flynn, ibid. 35, 1157 (1962). 

[6] R. Simha and L. A. Wall, J. Phys. Chem. 56, 707- (1952). 

[7] L. A. Wall, S. Straus, J. H. Flynn, D. Mclntyre, and R. Simha, J. Phys. Chem. 70, 53 (1966). 

[8] A. W. Smith and S. Aranoff, J. Phys. Chem. 62, 684 (1958). 

[9] P. Vallet, Compt. Rendus hebd. Seances Acad. Sci. (Paris) 249, 823 (1959). 
[10] C. D. Doyle, J. Appl. Polymer Sci. 5, 285 (1961). 
[11] J. H. Flynn and L. A. Wall, Polymer Letters 4, 323 (1966). 
[12] E. Jahnke, F. Loesch, and F. Emde, Tables of Higher Functions, 6th ed., p. 18 (McGraw-Hill Book Co. New York, 

N.Y., 1960). 
[13] O. Schlomilch, Vorlesungen iiber hoheren Analysis, 2d ed., p. 266 (Bruenschweig, 1874). 
[14] C. D. Doyle, J. Appl. Polymer Sci. 6, 639 (1962). 
[15] C. D. Doyle, Nature 207, 290 (1965). 
[16] R. W. Farmer, Wright-Patterson Air Force Base, Ohio, Aeronautical Systems Division Tech. Doc. Rept. ASD-TDR- 

62-1043, Part I (1963). 
[17] D. W. van Krevelen, C. van Heerden, and F. J. Huntjens, Fuel 30, 253 (1951). 
[18] J. P. Baur, D. W. Bridges, and W. M. Frassell, J. Electro-chem. Soc. 102, 490 (1955). 
[19] A. W. Coats and J. P. Redfern, Nature 201, 68 (1964). 
[20] A. W. Coats and J. P. Redfern, Polymer Letters 3, 917 (1965). 
[21] H. H. Horowitz and G. Metzger, Analytical Chemistry 35, 1464 (1963). 
[22] J. R. Soulen, Analytical Chemistry 34, 136 (1962)! 

[23] F. L. McCrackin, Polymer Physics Section, National Bureau of Standards, private communication. 
[24] L. Reich and D. W. Levi, Makromol. Chem. 66, lf)2 (1963). 
[25] C. D. Doyle, Makromol. Chem. 80, 220 (1964). 
[26] L. Reich, Polymer Letters 3, 231 (1965). 
[27] L. Reich, Polymer Letters 2, 621 (1964). 

[27a] T. Ozawa, Bui Chem. Soc. Japan 38 [11], 1881-1886 (Nov. 1965). 
[28] R. C. Turner and M. Schnitzer, Soil Science 93, 225 (1962). 
[29] R. C. Turner, I. Hoffman, and D. Chen. Canadian J. of Chem. 41, 243 (1963). 
[30] B. Kaesche-Krischer and H. J. Heinrich, Chemie-Ing.-Techn. 32, 598 (1960). 
[31] P. Murry and J. White, Trans. Brit. Ceram. Soc. 54, 204 (1955). 

522 



[32] H. E. Kissinger, Analytical Chem. 29, 1702 (1957). 

[33] H. E. Kissinger, J. Res. NBS 57, 217 (1956). 

[34] R. M. Fuoss, I. 0. Salyer, and H. S. Wilson, J. Polymer Sci. 2A, 3147 (1964). 

[35] P.Kofstad, Nature 179, 1362 (1957). 

[36] R. M. Barrer and G. C. Bratt, J. Phys Chem. Solids 12, 130 (1959). 

[37] A. E. Newkirk, Analytical Chem. 32, 1558 (I960). 

[38] D. A. Anderson and E. S. Freeman, J. Polymer Sci. 54, 253 (1961). 

[39] K. Traore, Soc. Chim. France, 5 serie, Memoirs, p. 274 (1964). 

[40] J. A. Magnuson, Analytical Chem. 36, 1807 (1964). 

[41] T. R. Ingraham and P. Marier, Canadian J. Chem. Eng. 42, 161 (1964). 

[42] H. C. Anderson, J. Polymer Sci. 6C, 175 (1965). 

[43] H. L. Friedman, ibid., 6C, 183 (1965). 

[44] P. K. Chatterjee, J. Polymer Sci. 3A, 4253 (1965). 

[45] J. H. Flynn, Chem. Eng. 69, 137 (1962). 

[46] E. S. Freeman and B. Carroll, J. Phys. Chem. 62, 394 (1958). 

[47] P. Lumme, Suomen Kemistilehti B32, 198 (1959). 

[48] P. Lumme and J. Peltonen, ibid. B37, 162 (1964). 

[49] H. C. Anderson, Makromol. Chem. 51, 233 (1962). 

[50] S. Igarashi and H. Kambe, Bui. Chem. Soc. Japan 37, 176 (1964). 

[51] H.-J. Frentz, Elektrotechnische Z. A78, 156 (1957). 

[52] L. C. Whitman and W. W. Whitman, Amer. Inst. Elect. Eng. Trans. 78, Pt. I, 308 (1959). 

[53] L. Reich, H. T. Lee, and D. W. Levi, J. Applied Polymer Sci. 9, 351 (1965). 

[54] H. C. Anderson, Polymer Letters 2, 115 (1964). 



(Paper 70A6-422) 



523 



